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ABSTRACT 

Centroid positions of peaks identified in weak lensing mass maps often show offsets 
with respect to other means of identifying halo centres, hke position of the brightest 
cluster galaxy or X-ray emission centroid. Here we study the effect of projected large- 
scale structure (LSS), smoothing of mass maps, and shape noise on the weak lensing 
peak positions. Additionally we compare the offsets in mass maps to those found in 
parametric model fits. Using ray-tracing simulations through the Millennium Run 
A'^-body simulation, we find that projected LSS does not alter the weak-lensing peak 
position within the limits of our simulations' spatial resolution, which exceeds the 
typical resolution of weak lensing maps. We conclude that projected LSS, although a 
major contaminant for weak-lensing mass estimates, is not a source of confusion for 
identifying halo centres. The typically reported offsets in the literature are caused by 
a combination of shape noise and smoothing alone. This is true for centroid positions 
derived both from mass maps and model fits. 
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1 INTRODUCTION 

Weak gravitational lensing is one of the most important 
tools to study clusters of galaxies. Its sensitivity to both 
luminous and dark matter is a distinct advantage over 
other methods. Observations of a clusters' tidal shear 
field imprinted on the observed ellipticities of faint back- 
ground galaxies allow for "parameter-fre e " mass reconstruc- 
tions (see e.g. iKaiser fc Squired 1 19 93'; 'Soitz fc S chneidej 
200ll for methods, and IClowe fcSchneid cr 200ll. l2002f 

2010l 



Dietrich etHI |2009| ; llsrael et al 



for applications). 
Such two-dimensional maps often show offsets between the 
position of the observed peaks in surface mass density and 
other indications of a galaxy cluster's centre, like the posi- 
tion of the brightest cluster galaxy or the centroid of the 



X-ray emission (e.g., Dietrich et al.l |2005| ; IWittman et al.l 
l2006l : IOkabe et al.ll201o ir Such offsets can be of either astro- 
physical origin or caused by random or systematic noise. An 
example of the latter could be missing data due to bright 
stars or other masked image defects. 

Two examples of clusters where the observed sep- 
aration between weak-lensing and X-ray centroid has 
been attr ibuted to astrophys ical causes are the "Bullet 
Cluster" (|Clowe et all l2006al ) and MACS J0025.4-1222 
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IIBradac et al ] |2008l '). These are systems in which high-speed 
collisions separated the coUisionless dark-matter and galaxy 
components from the coUisional intra-cluster medium. 

In other cases the spatial discrepancy has been de- 
termined to be consistent with the centroid shifts in- 
duced by the shape noise of the background galaxies (e.g., 
IClowe et al.l l2006bl ) . It has also been noted that the ob- 
served centroid offsets depend on the halo concentration, 
with more centrally concentrated halos showing smaller off- 
sets in convergence maps (?). 

Weak lensing is susceptible to the entire mass along 
the line-of-sight and it is now well established that pro- 
jected large-scale structure (LSS) and halo triaxiality 
are significant contributions to the total error budget in 
weak - lensin g m ass estimates of galaxy clusters jHoekstral 



ik-lensmg m ass estimates ot galaxy clusters lllloekstral 
H. l2003l; ?; iDodelsonl |2004| ; iBecker fc Kravtsovl l20ia 
Corless fc Kind 2007 ). In addition to infiuencing weak- 
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lensing mass-estimates, irregular halo shapes and projected 
LSS could potentially shift the centroids of the observed 
surface mass density distribution away from the true halo 
centre. We will summarily refer to t hese effects as pr o jected 
LSS throughout this work. E.g., iLombardi et al.l j2005l ) 
speculated that the projection of groups close to the line- 
of-sight might be the cause for the 8" or 94/!,^^kpc offset 
they observed between the weak lensing and optical/X-ray 
centroids of RCDS 1259.9 — 2927. Our aim in this paper is 
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to study the importance of this effect and its relative im- 
portance to centroid shifts induced by shape noise and the 
inevitable smoothing in mass reconstructions. Knowledge of 
all these effects is required to make robust statements about 
the nature and astrophysical significance of observed peak 
offsets. 

The outline of the paper is as follows. We describe the 
simulations and how we identified weak lensing peaks in 
them in Sect. [2] We present the results of matching halo 
centres to lensing peak positions in Sect.O Finally, we sum- 
marise and discuss our findings i n Sect.lH ^ye use standard 
lensing notation throughout (e.g. ISchneidedlioOfil) . 



2 METHODS 

To study the spatial offsets between halo centre and lensing 
peak position, we simulated weak-lensing observations of 
galaxy clusters using A^'-body and ra y-tra cing simul ations . 
We used the Millennium Run (MR, ISpringel etlill boosl ) 
simulation. The MR is a large dark matter A'^-body simula- 
tion with 2160^ particles in a 500/i~^Mpc box simulating 
structure formation from z — 127 to the present epoch in 
a flat ACDM Universe. The assumed cosmological param- 
eters are: a matter density Sim ~ 0.25, a cosmological con- 
stant with energy density Q,a ~ 0.75, a Hubble constant of 
h — 0.73 (in units of 100kms~^Mpc~^), and an initial den- 
sity power spectrum normalisation erg = 0.9. The Plummer 
equivalent smoothing scale in the MR is 5 kpc. 

Cluster halos in the MR were identified by a friend- 
of-friends algorithm. Cluster centres were identified with 
the position of the minimum of the gravitational potential. 
Cluster masses were obtained by measuring the mass inside 
spheres around the cluster centre with mean density 200 

times the critical cosmic density. 

The ray-tracing algorithm described in iHilbert et al.l 

(l2009l ) was used to calculate the apparent positions of the 
cluster centres and the convergence in 512 4 x 4deg^ fields 
of view. In these simulations, the matter distribution is cre- 
ated by periodic continuations of the simulation boxes at 
increasing redshifts to create backward light cones. Matter 
is divided into slices of thickness ~ 100 h^^ Mpc perpendic- 
ular to the line-of-sight and angled with respect to the box 
sides to avoid the repetition of structures along the line- 
of-sight. Light rays are shot from redshift zero to 2: = 3.06 
through these lens planes. The solid angle subtended by 
these simulations is small enough to neglect the effects of 
sky curvature. 

The products of these ray-tracing simulations are noise- 
free maps of the dimensionless surface mass-density k and 
complex shear 7. This noise-free case allows us to study 
the impact of LSS projections alone, without being affected 
by smoothing and shape noise, which are always present 
in mass reconstructions derived from astronomical observa- 
tions. 

2.1 Mass maps 

Already iKaiser fc Squired (|l993h noted that a mass recon- 
struction obtained from an unfiltered shear field would 
have infinite variance. Smoothing the shear field with a 
low-pass filter suppresses the noise but also amounts to a 
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Figure 1. Physical offsets between 512 halo and lensing peak 
positions for a single background redshift population at z = 1.08. 

smoothing of the reconstructed k field. F or a Gaussian fil- 
ter PC exp{—6^/2a^) the variance in k is (|Kaiser fc Squirei 
I1993I ) 

where is the intrinsic 2-d ellipticity dispersion and n the 
number density of background galaxies. This smoothing also 
leads to correlated noise in the mass reconstruction. The 
noise power spectrum is given by (|Lombardi fc Bertinlll99^ ) 

The noise on large fields is rotation and translation invari- 
ant, and therefore the covariance between modes at different 
k is zero. Maps of correlated noise were created by draw- 
ing independent complex Gaussian random variables with 
variance given by eq. ((l| and applying an inverse Fourier 
transform to the arrays created in this way. 

In this work we used three values for the number 
density of background galaxies n — {10, 30, 80} arcmin"^ 
and two smoothing scales — {45", 90"}. The two 
smaller number densities correspond roughly to typical 
groun d-based observations with 2 m to 8 m cla, ss telescopes 
(e.g., iDietrich et aD l2007l : lOkabe et all l2010h . while the 
higher number de nsity represents space-based data (e.g., 
IClowe et al.ll2006al l . The large smoothing scale was used for 
the lowest number density, while the small smoothing scale 
was applied in the case of the two larger number densities. 

Peaks in K-maps were detected with a connected- 
component labelling of pixels above a detection threshold. 
We used the 8-connectivity, i.e., we consider all pixels that 
are connected via the sides, edges, or corners of a square as 
one structure. Additionally, we imposed the condition that 
each peak must have at least 3 pixels above the detection 
threshold. The pixel with the highest value is considered to 
be the peak location. We compar ed this peak location to 
the one derived by SExtractor's (jBertin fc Arnorit3ll996l ') 
fit of a bivariate Gaussian. The two methods give results 
that generally agree to within a fraction of pixel. This de- 
viation is too small to be a concern for us. Halos were then 
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matched to their nearest neighbour in the hst of detected 
peaks. All halos with a mass of more than 10^* Mq and 
2 < 1.1 were considered in this matching. 

2.2 Shear catalogues and model fits 

Often weak-lensing studies of clusters include paramet- 
ric model fits to the tangential s hear profiles (e.g., 
iDietrich et "alll2009l : [Okabe et al.ll2010l l. in which the cen- 
troid position ca n be treated as free parameters (e.g., 
lOguri et al.ll2010h . Since no smoothing is applied to the 
shear field, shape noise and LSS projections are the only 
source of noise. To compare with the centroid positions 
obtained from mass maps, we created mock catalogues of 
galaxies with number densities n = {10, 30, 80} arcmin"^ 
and intrinsic ellipticities drawn from a Gaussian with 2- 
dimensional ellipticity dispersion = 0.38, truncated at 
|e| — 1. The reduced shear g, computed from the ray- 
tracing, was applied to these intrinsic ellipticities using the 
standard relation (?) 



e*( = )-|-9- 



(3) 



where the asterisk denotes complex conjugation, and the 
superscript (s) denotes intrinsic source properties. 

We fitted NFW (?) models to the shear estimated from 
the galaxy catalogue within a radius of 15' around the halo 
centre positions. We fixed the NFW concentration parame- 
ter to follow the ? mass-concentration relation, leaving three 
free parameter, the halo centre {xc,yc) and the halo mass 
771200- Fits were performed using the maximum likelihood 
estimator of ?. 



3 RESULTS 

3.1 Noise-free case 

Figure [T] shows the apparent and physical offsets between 
centres of halos and associated peak positions in noise free 
K-maps. While the matching between halos and lensing 
peaks was performed for all cluster sized halos at redshifts 
below 1.1 as described in the previous section, we made a 
redshift dependent mass cut in the selection of halos that 
enter our comparison. We selected a fiducial minimum clus- 
ter mass of Mfld — 2 x 10^"* Mq at redshift z^d = 0.3. 
The minimum mass is governed by the lensing efficiency and 
scales with the halo redshift z as 

X>(0, Zfld)©(2fld, Zhgk) 



Mmin{z) = A/fld- 



©(0,2)0(2, 2bgk) 



(4) 



where X>(-, ■) is the angular diameter distance between two 
redshifts, and 2bgk is the redshift of the background source 
population. Ifere the background sources were assumed to 
be at a single redshift 2bgk ~ 1.08, which is the redshift of 
the first lens plane at 2 > 1. 

Obviously, almost all lensing peaks are projected to the 
immediate vicinity of the halo centre. Virtually no associa- 
tions are made on scales > 15 kpc, which are beyond the 
resolution limit of the Millennium Run simulation. Visual 
inspection of the few large-distance matches confirmed that 
their large offsets are not caused by LSS projections. They 




rf/arciniii 

Figure 2. Apparent offsets in smoothed K-maps. The red/biue 
(rising/faffing hatches) histograms are for smoothing scales of 
45" and 90", respectively. The thin solid green histogram is the 
unsmoothed case from Fig. [T] shown here for comparison. 



are undetected halos erroneously matched to a nearby lens- 
ing peak or their surface mass density signal was merged 
into a double peaked structure when it overlapped with 
a more significant foreground or background object. Our 
peak finder does not attempt to deblend such objects into 
multiple peaks, although a human observer would correctly 
identify the smaller peak as caused by the halo. Thus, we 
do not consider these cases to be LSS projections for the 
purpose of this work. Furthermore, in optical multi-colour 
data, both clusters would be readily identifiable as separate 
entities. 

While a source redshift of 2 ~ 1 is typical for the mean 
redshift of targeted ground-based lensing observations, the 
tail of the redshift distribution extends to higher values, and 
space based observations have the potential to probe LSS 
at significantly higher redshifts. We therefore matched the 
same set of halos to convergence peaks detected in maps 
with a single source plane at z — 3.06, the highest redshift 
in our ray-tracing simulations, in an effort to maximise the 
impact of LSS projections. 

The result is virtually unchanged. The 90% percentile 
increases from 2.0 kpc to 5.6/i~^kpc and a few more 
halos have barely resolved offsets, but still 478 out of the 
512 halos have offsets of less than 15/i~^kpc, down from 
503 in the case of the lower source redshift. 



3.2 Smootliing and shape noise 

We now turn to the influences of smoothing and shape noise. 
The effect of smoothing is shown in Fig. [2] 474 (434) halos 
fulfilling our selection criterion in eq. Q could be matched 
to peaks detected in K-maps smoothed with Gaussians of 
width 45" (90") within a distance of 5'. Again, the vast 
majority of large separation pairs in the tail of the distri- 
bution is due to erroneous associations. Ignoring outliers 
at d > 2', we find that almost 90% of all associations are 
made within one half of the smoothing scale. The observed 
offsets are, however, significantly larger than in the absence 
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Figure 3. Apparent offsets in smoothed noisy K-maps. The three 
panels are for different smoothing scales and number densities as 
denoted in their legends. 



Figure 4. Apparent centroid offsets in NFW model fits for three 
different number densities. 



of smoothing, as is shown by the comparison with the un- 
smoothed case in Fig. (2] Thus, smoothing alone, without 
the added influence of shape noise, causes shifts of centroid 
positions with respect to the halo centres. 

We summarise properties of the offset distributions in 
Table [T] Here the mode of the distributions is estimated 
from the maximum of a Gaussian kernel density estimate 
of the offset distribution. It is worth pointing out that for 
one-sided distributions whose mode is very close to zero, this 
method has a systematic bias towards higher values. Never- 
theless, the values we flnd are generally in good agreement 
with naive mode estimates obtained from looking at the his- 
togram plots, with the exception of the LSS only case. All 
other quantities are computed directly from the measured 
offsets. 

Finally, Fig. [3] shows the observed offsets in smoothed 
maps with shape noise added. As described in Sect. [2] we 
assumed combinations of number densities and smoothing 
scales of (n/arcmin-^ cts/" = {(10, 90), (30, 45), (80, 45)}. 
As one would expect, the apparent offsets are larger than 
in the noise free case. The distribution's mode roughly 
triples for n = 10 arcmin"^ and approximately doubles for 
n = {30, 80} arcmin"^. The distributions also have notice- 
ably longer tails than in the noise free case. Again, the dis- 
tribution properties are summarised in Table [T] 



3.3 NFW model fits 

The centroid offset distributions obtained from fitting NFW 
models to the shear catalogues, as described in section \2l2l 
are shown in Fig. |4] and the distribution properties are sum- 
marised in Table [T] The offsets are smaller than those ob- 
served in K-maps for the low number density case only. For 
n = {30, 80} arcmin"'^ the offsets are comparable, if not 
slightly worse than in noisy mass maps. 



4 SUMMARY AND DISCUSSION 

We used ray-tracing simulations through the Millennium 
Run A^-body simulation to study the impact of projected 
large-scale structure on the observed peak positions of ha- 
los in weak-lensing maps. The main result of this work is 
that such projections do not lead to shifts of the maximum 
surface mass density away from the halo centre. This is true 
even for the extreme case of a single source redshift at z = 3. 
Thus, this is a robust result applicable to both ground- and 
deep space-based lensing observations. By making the mass 
selection in eq. Q dependent only on the lensing efficiency 
in our mock survey, we neglected selection biases that will 
effect real galaxy cluster catalogues. Real clusters are often 
studied in detail precisely because they reveal striking lens- 
ing features, offsets between dark and luminous matter, or 
other unusual features. Furthermore, mergers, projections, 
and alignment of triaxial halos with the line-of-sight will 
inevitably boost low mass halos to appear to be above the 
mass selection. While this could in principle be a concern 
for the applicability of our results to a highly biased cluster 
sample, the virtually complete absence of significant offsets 
caused by LSS projections across our entire sample suggests 
that this is not a concern. 

The negligible effect that LSS projections have is ex- 
emplified by the case of RCDS 1259.9 - 2927 at z = 1.24, 
which we mentioned in the introduction. Its lensing cen- 
troid has an offset of 8". This corresponds to a projected 
distance of 48/i~^kpc and is significantly larger than any 
of the offsets seen in Fig. [1] Thus, this offset is most likely 
not caused by projected LSS. We emphasise that, although 
Table [1] suggests that an offset of 8" is consistent with LSS 
projections at the ~ 95% confidence level, we see such large 
offsets only for lower redshift halos. The relevant compari- 
son here is with the physical offsets reported in Fig. [T] This 
is supported by looking at the offsets caused by LSS as a 
function of redshift. We grouped the 512 halos into four 
equally large groups of increasing redshifts and report the 
median and 95th percentile of their respective offset distri- 
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Table 1. Properties of the offsets distributions in different K-maps and in the model fits as a function of number density n and smoothing 
scale cTs. 
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Table 2. Offsets caused by LSS projections as a function of halo 
rodshift. 
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butions in Table [21 The median shows a slight increase with 
redshift but still well below the resolution of our simula- 
tions. The tails of the distributions show no clear trend. We 
conclude that the offsets caused by LSS are limited in phys- 
ical scales to those scales that are much smaller than the 
redshift independent smoothing scale for all halos studied 
in this work. 

Smoothing an asymmetric mass distribution with a 
symmetric kernel naturally leads to centroid shifts. We 
showed that these occur on scales less or equal to one half of 
the smoothing scale. Adding shape noise leads to additional, 
significantly larger offsets. 

For typical shallow ground-based observations with 
number density lOarcmin^'^ and smoothing scale 90", the 
mode of the offset distribution is at ~ 44" and the mean 
of the distribution is at 52". These values are compati- 
ble with tho se found by blind searches for galaxy clusters 
(iGavazzi fc So ucail 2007; Schirmcr et al. 2007; Gc Uer et al] 
I2OIOI ). All these values were obtained by taking only those 
halos into account which could be matched to a weak- 
lensing peak within 5'. Including the long tail of matches at 
larger separations of course leads to the inclusion of mostly 
spurious matches. The 75% percent ile of the full distribut ion 
is at 2'.3, reproducing the result of iDietrich et al.l (|2007l ). 

Obviously the offsets become smaller if a higher number 
density of background galaxies suppresses the shape noise. 
However, even with the depths achieved by current space- 
based data sets, the offsets are significantly larger than in 
the absence of shape noise. E.g., for n — 80arcmin~^ the 
mode of the offset distribution is at 15", almost twice as 



large as in the noise- free case where it is at 8". The situation 
is slightly better for the tail of the distribution where the 
95th percentile is still larger by a factor of 1.5. 



Our findings show that the significance of peak offsets 
can often be reliably estimated by bo otstrapping f rom th e 
shear catalogue, as it was e.g. done bv lClowe et al.l(l2006al ). 
Such a bootstrapping procedure does not include the sys- 
tematic centroid shifts due to smoothing but this effect is 
often smaller than the offsets typically caused by shape 
noise. The resulting underestimation of offset distances can 
be rectified by taking into account that smoothing leads 
to additional offsets of up to half the smoothing scale at 
90% confidence. Applying this to the example of the Bullet 
Cluster, we conclude that the contribution of smoothing is 
absolutely negligible and the bootstrapping procedure pro- 
vided a reliable estimate of the offsets' significances. 



lOguri elal] l|2010l l studied the offsets of the weak lens- 
ing signal of 25 massive clusters in the LoCuSS sample from 
their brightest cluster galaxies (BCG). When comparing 
their results to our NFW model fits, it is important to re- 
member that they have no knowledge of the true halo centre 
as we have in the MR but must instead rely on the BCG 
position as an estimate for the halo centre. This is subject 
to possible misidentifications of the BCG or displacements 
of the BCG with respect to the dark matter h alo centre, 
as discussed in detail in [Johnston et al] (|2007l '). However, 
lOguri et al.l conclude that their result is consistent with a 
10% fraction of clusters wi th BCG m isidenitifications, in 
agreement with the work of [Johnston et al. (2007;). The re- 
maining offsets are all smaller than 1'. This is comparable 
to what we find in the middle panel of Fig. (4] which has ap- 
proximately the same number density as the LoCuSS data 
set. However, our results also suggest that ~ 10% of all off- 
sets are larger than 1'. This could very well masquerade as 
BCG misidentifications, which should occur with the same 
frequency. 
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